Aproximativ 20% de teste la metode de calcul se refera la rezolvarea sistemelor de ecuatii liniare si matrici Tema respectiva se studiaza în detalii în cursul matematicii în clasele a XI-a - a XII-a de liceu Din aceste considerente vom aminti numai notiunile de baza Se considera un sistem de m ecuatii liniare cu n necunoscute x1, x2, , xn, care poate fi scris astfel: (1) Aceasta scriere presupune ca ecuatiile sistemului sunt numerotate cu numerele 1, 2, , m Coeficientii pe lînga necunoscute au cîte 2 indici, primul indica numarul ecuatiei, iar al doilea- numarul necunoscutei De exemplu, aij este coeficientul pe lînga necunoscuta xj în ecuatia i Vom scrie coeficientii pe lînga necunoscutele sistemului (1) în forma de matrice: Asa matrice poarta denumirea de matrice a coeficientilor sistemului Daca vom completa matricea cu o coloana alcatuita din termenii liberi al sistemului, vom capata o matrice noua, care poarta denumirea de matrice extinsa a sistemului: Numerele ?1, ?2, , ?n se numesc solutie a sistemului (1), daca substituirea necunoscutelor x1, x2, , xn respectiv prin numerele?1, ?2, , ?n transforma toate ecuatiile sistemului în egalitati Sistemul se numeste compatibil daca exista cel putin o solutie a acestuia În caz contrar sistemul este incompatibil Un sistem compatibil poate avea o singura solutie (sistem determinat) si multe solutii (sistem nedeterminat) Vom descrie în continuare metoda Gauss de rezolvare a sistemului de n ecuatii liniare cu n necunoscute Calculul solutiei sistemului prin metoda Gauss consta în eliminarea succesiva a necunoscutelor prin transformari echivalente a sistemului la o forma cu matrice superioara triunghiulara (etapa mersului direct) Necunoscutele se calculeaza la etapa mersului invers Vom pastra matricea extinsa a sistemului în memoria calculatorului în forma unui tablou bidimensional de numere reale, Asadar, termenii liberi bi se noteaza în continuare ai,n+1 Etapa mersului direct consta din n-1 pasi Pasul 1 Ca scop se pune eliminarea necunoscutei x1 din ecuatiile cu numerele i = 2, , n Presupunem ca coeficientul a11 ? 0 În caz contrar trebuie sa alegem din ecuatiile plasate mai jos acea ecuatie la care ai1 ? 0, si sa schimbam cu locurile prima ecuatie si ecuatia i La toti pasii care urmeaza aceasta operatie se repeta, în caz de necesitate (Numarul de schimbari efectuate este important, pentru ca fiecare schimbare de ecuatii influienteaza semnul determinantului matricei ) Vom calcula ci1 = ai1/a11 (i = 2, 3, , n) si vom scadea consecutiv din ecuatiile 2, 3, , n prima ecuatie înmultita respectiv la c21, c31, , cn1 Aceste operatii transforma în zero coeficientii pe lînga x1 în toate ecuatiile cu exceptia primei Ca rezultat capatam un siste echivalent de ecuatii a11x1 + a12x2 + a13x3 + + a1nxn = b1 , a22(1)x2 + a23(1)x3 + + a2n(1)xn = b2(1) , a32(1)x2 + a33(1)x3 + + a3n(1)xn = b3(1) , an2(1)x2 + an3(1)x3 + + ann(1)xn = bn(1) în care aij(1) ? bij(1) se calculeaza asa: aij(1) = aij - ci1a1j , bi(1) = bi - ci1b1 La fel procedam si cu celelalte necunoscute Vom descrie pasul k Pasul k În conditii de akk(k-1) ? 0 calculam cik = aik(k-1) / akk(k-1) (i = k + 1, , n) si scadem din ecuatiile (k + 1), , n ecuatia k obtinuta la precedentul pas si inmultita respectiv la ck+1,k, ck+2,k, , cnk Dupa pasul (n - 1) de eliminari obtinem urmatorul sistem de ecuatii: a11x1 + a12x2 + a13x3 + + a1nxn = b1 , a22(1)x2 + a23(1)x3 + + a2n(1)xn = b2(1) , a33(2)x3 + + a3n(2)xn = b3(2) , ann(n-1)xn = bn(n-1) cu matricea A(n-1) de coeficienti superior triunghiulara Aici calculele etapei mersului direct se încheie Etapa mersului invers Din ultima ecuatie a sistemului calculam xn Substituim xn prin valoarea gasita în toate ecuatiile începînd cu a (n-1) pîna la prima calculam consecutiv xn-1, xn-2, , x1 Necunoscutele se calculeaza dupa urmatoarele formule: xn = bn(n-1) / ann(n-1), xk = (bn(k-1) - ak,k+1(k-1)xk+1 - - akn(k-1)xn) / akk(k-1), (k = n - 1, , 1) Propunem mai jos programul Pascal, care realizeaza metoda Gauss Program Gauss; Type Matricea = array of real; vector = array of real; Var n : integer; a : matricea; x : vector; G : Boolean; n - numarul de ecuatii a - matricea coeficientilor sistemului x -vectorul necunoscutelor g - parametrul, care semnalizeaza faptul, este compatibil sau nu sistemul de ecuatii Procedure Citire ( var a : matricea; var n : integer); Var i, j : integer; Begin readln(n); for i :=1 to n do for j :=1 to n+1 do read(a[i,j]); end; Introducerea numarului de ecuatii n si a coeficientilor sistemului ai,j de la tastatura Procedure Direct ( var A : matricea; N : integer; var G : boolean); Var p,q,r, i, j : integer; c : real; Procedure Schimb (var a :matricea; n, p: integer; var G:boolean); Var i, j : integer; t : real; begin G:=true; i :=p+1; while (a[i,p]=0) and (i 0 then begin C:= a[I,1]/a[j,1]; { C - coeficientul de proportionalitate } For k:=2 to n+1 do {acum controlam celelalte perechi din liniile i si j} If a[j,k]<>0 then if a[i,k]/a[j,k]<>C then Pr:=False end else Pr:=false; if Pr then Goto 1; {Întrerupem procesul de control daca gasim doua linii proportionale } end; 1: Prop:=pr; end; Metoda lui Cramer Se considera un sistem de n ecuatii liniare cu n necunoscute x1, x2, , xn, care poate fi scris astfel: (1) Fie, determinantul matricei sistemului este nenul Vom nota: ; ; Sistemul (1) se poate rescrie sub forma unei ecuatii matriciale: AX=B Fie d=det A determinantul sistemului si dj, 1 ? j ? n, determinantul matricei Aj care se obtine din A prin înlocuirea coloanei j prin coloana B De exemplu: , s a m d Regula lui Cramer: Cu notatiile de mai înainte, daca d=det A este nenul (cu alte cuvinte, sistemul este compatibil determinat), atunci sistemul (1) are o solutie unica, anume: Aceste formule poarta denumirea formule lui Cramer program Cramer; uses crt; const nmax=20; type matr=array of real; vect= array of real; var a:matr; b,x :vect; d :real; f :boolean; i,n : integer; a - matricea coeficietilor sistemului; b - vectorul termenilor liberi; x - vectorul necunoscutelor; d - determinantul matricei coeficientilor; f - variabila booleana care semnalizeaza existenta solutiei; n - numarul de ecuatii în sistem procedure cit(var n:integer;var a:matr;var b:vect); var i,j:integer; begin Writeln('Tastati n numarul de ecuatii=>'); readln(n); writeln('Tastati el matricii coeficientilor=>'); for i:=1 to n do for j:=1 to n do begin write('a[', i , ',' , j ,']='); read(a[i,j]); end; writeln('Tastati terminii liberi=>'); for i:=1 to n do begin write('b[', i ,']='); read(b[i]); end; end; Procedura cit realizeza citirea datelor initiale: dimensiuninea matricei coeficientilor, coeficientii si termeni liberi ale sistemului de ecuatii function det(a:matr):real; var i:integer; aux,p:real; Calculul determinantului matricei coeficintilor sistemului de ecuatii procedure triungh( var f:boolean); var i,j,k,i1:integer; p:real; Procedura triungh modifica matricea coeficientilor la o forma superior triunghiulara prin transformari echivalente procedure permut(k,l:integer); var i:integer; c:real; begin for i:=k to n do begin c:=a[k,i]; a[k,i]:=a[l,i]; a[l,i]:=c; end; end; Procedura Permut inerschimba cu locul doua linii în matricea coeficientilor în cazul cînd elementul diagonal în linia i este egal cu 0 begin f:=true; for i:=1 to n-1 do begin i1:=i; while (a[i1,i]=0) or (i1>n) do i1:=i1+1; if i1>n then f:=false else begin if i1<>i then begin permut(i,i1);aux:=-aux;end; for j:=i+1 to n do begin p:=-a[j,i]/a[i,i]; for k:=i to n+1 do a[j,k]:=a[j,k]+a[i,k]*p; end; end; end; end; Fiecare schimbare de linii cu locurile modifica semnul determinantului Variabila intermediara aux retine aceasta schimbare Begin aux:=1; triungh(f); if f then begin for i:=1 to n do aux:=aux*a[i,i]; det:=aux; end else det:=0; end; Valoarea TRUE a variabilei f arata ca sistemul este compatibil determinat si are o singura solutie, în caz contrar determinantul matricei coeficientilor este egal cu 0 Determinantul matricei se calculeaza ca produsul elementelor de pe diagonala principala, datorita faptului ca toate elementele sub diagonala principala sînt egale cu 0 si, prin urmare, ceilalti termeni ai determinantului sînt egali cu 0 procedure kramer(a:matr;var x:vect); var i:integer; d,d1:real; var aa:matr; Procedura Kramer realizeaza regula lui Cramer procedure subst(var a:matr;b:vect; l:integer); var i:integer; begin for i:=1 to n do a[i,l]:=b[i]; end; Modificarea matricei coeficientilor prin substituirea coloanei l prin vectorul B Begin d:=det(a); for i:=1 to n do begin aa:=a; subst(aa,b,i); d1:=det(aa); x[i]:= d1/d; end; end; aa - matricea intermediara, în care toate elementele, cu exceptia celor din coloana i, coincid cu elementele matricei coeficientilor d1 - determinantul matricei aa, obtinute din a prin substituirea coloanei i prin vectorul B x - vectorul necunoscutelor, valorile carora se calculeaza dupa regula lui Cramer : di/d Begin clrscr; cit(n,a,b); d:=det(a); writeln('det=',d:6:2); if d=0 then writeln('Sistemul de ecuatii nu are solutii') else begin kramer(a,x); writeln('solutia prin metoda Kramer:'); for i:=1 to n do writeln(x[i]:6:2); end; end Partea operativa a programului principal Aici initial apelam procedura cit pentu a introduce dimensiunea si coeficientii sistemului de ecuatii Daca determinantul sistemului este egal cu 0, atunci sistemul nu are solutii, altfel procedura Kramer ne returneaza vectorul X - solutia sistemului Afisam la ecran componentele vectorului X 